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In the presence of strong random scattering the behavior of particles with degenerate spectra is 
quite different from Anderson localization of particles in a single band: it creates geometric states 
rather than confining the particles to an area of the size of the localization length. These states are 
subject to a Fokker-Planck dynamics with universal drift velocity and disorder dependent diffusion 
coefficient. This behavior has some similarity with the unidirectionally propagating edge states in 
quantum Hall systems. 
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Random scattering of wave-like states (electrons, 
photons or acoustic waves) leads either to diffu¬ 
sion for weak random scattering or to Anderson 
localization (AL) in the presence of strong random 
scattering. AL is a phenomenon where diffusion 
or propagation is suppressed because random scat¬ 
tering confines the modes to a finite region, whose 
size is characterized by the localization length [II- 
. The characteristics of AL can be observed in a 
pure form for photons [3-Q due to the absence of 
an additional particle-particle interaction. 


index and which transforms H into —H* as 

H = -UH*U^. (1) 

In this paper we will show that ray modes can 
be created spontaneously by strong random scat¬ 
tering in systems with spectral degeneracies based 
on the relation (HD. Starting point is the transition 
probability for a particle, governed by the random 
Hamiltonian H, to move from the site r' on a lat¬ 
tice to another lattice site r within the time t: 


Another interesting phenomenon in disordered 
systems is the quantum Hall effect, which is char¬ 
acterized by plateaux in the Hall conductivity 
0- The latter have been attributed to propagat¬ 
ing edge modes in the otherwise Anderson local¬ 
ized bulk system Q. Edge modes also exist in 
gapped systems without any random scattering. 
A simple case is a massive 2D Dirac Hamiltonian, 
whose mass sign changes by crossing an edge in 
y-direction: The mass is —m for a; < 0 and m for 
a; > 0. The gap 2|m| > 0 prevents the system to 
generate any other extended state than that along 
the edge. This mode decays exponentially when 
we go away from the edge in a:-direction. The ap- 
pearence of such states can be realized in a pho¬ 
tonic crystal with Faraday effect @ . 

It was recently observed that similar one¬ 
dimensional modes can also be generated sponta¬ 
neously by strong scatterin g in 2D systems without 
the existence of any edge |lCll . ini. These states, 
which will be called ray modes subsequently, are 
created in systems with a generalized particle-hole 
symmetry. The latter implies spectral degeneracies 
[ll|. This is caused by the fact that the Hamilto¬ 
nian H acts on states Abrj that depend on space 
coordinates r and on an additional spinor index j, 
typically with values j = 1,2. Then there exists 
a unitary matrix U that acts only on the spinor 


where (r, j\...\r', j') is the quantum average and 
is the average with respect to randomly dis¬ 
tributed disorder, r and r' refer to real space coor¬ 
dinates and the indices j, j' refer to different bands 
of the system. Prr'ir) is a fundamental quantity 
from which we can obtain transport and localiza¬ 
tion properties 

In the following we will focus on 2D Dirac and 
3D Weyl particles and on particles on the square 
lattice with tt flux. For all these models the Hamil¬ 
tonian reads in sublattice representation 

H = a-H + aoV , (3) 

where aj {j = 0,1,2,3) are Pauli matrices with 
the 2x2 unit matrix ag and a random potential 
V with mean zero and variance {VrVr')d = gSrr'- 
The part with cti, a 2 provides scattering between 
different values of j, which will be crucial for the 
subsequent discussion. 

The relation (HD is satisfied, for instance, for 
the block-diagonal matrix of the 3D gapless Weyl 
Hamiltonian diag{aipi + CT 2 P 2 + + Vao,o'ipi + 

CT 2 T 2 — CTspa — Hcto) with 




CTi 

0 



( 4 ) 
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where Pj is the momentum operator with p* = 
—Pj- Another example is the massive 2D Dirac 
Hamiltonian diag{aipi +a 2 P 2 + o' 3 m + Vao, aipi + 
(T 2 P 2 + o' 3 TO —Hcto), which also obeys condition ([1]). 

It was shown in Ref. [13 that the Fourier com¬ 
ponent Prr'{ie) of Prr'ir) agrees for large distances 
\r — r'\ with a correlation function of a random- 
phase model, defined by the random matrix 

Crr' — ‘^drr' ^ ^ hj-j j-'j' ^ ^ j/ ,^/fj// , (5) 

j,j' j" ,r" 

where the propagator 


there is an invariance of C with respect to a global 
phase change arj —> cxrj + (j)- This implies that the 
mean-field distribution depends only on the rela¬ 
tive phase difference A = di — d 2 - Its value is 
determined by the condition maxA log(|Cp(A)|), 

where C'p(A) are the Fourier components of C with 
uniform phases. This condition may have several 
solutions A; with the same maximum, such that 
we have to sum over all of them: 

det C « ^ exp 

i 


log(4(A0) 


( 10 ) 


Kr' = aoSrr' + 2ir]{'H - ifjaQ)~^, (6) 

depends on the random phase Hamiltonian 

nrj,r'f = . (7) 


It turns out that the Fourier components (7p(A/) 
can be written as 

C'p(A) = 2 — TTo — s • if (11) 


H = {H)ii is the average Hamiltonian and rj is the 
scattering rate while fj = p + e. p can be consid¬ 
ered as an empirical parameter or can be calculated 
self-consistently from the self-energy of the average 
one-particle Green’s function {{H — z)~^)d HI . In 
any case, it increases with variance g of the random 
potential. 

In the limit e —> 0 the propagator h is unitary 
since 


hh'^ = 1 — 4ep{'H'^ + p^ao) ^ . (8) 


Then there is the following asymptotic relation for 
large scales between the random phase model and 
the Fourier components of the average transition 
probability (l^ : 


Prr* (fe) 


{C-] det C)a 
{detC)a 


(9) 


where the brackets (...)a mean integration with re¬ 
spect to the angular variables {0 < arj < 27r}, nor¬ 
malized by 27r. These random angles represent the 
relevant part of the disorder fluctuations in terms 
of long-range correlations. In other words, there is 
a mapping from the original random Hamiltonian 
in Eq. ([3|) to the random phase Hamiltonian TL 
that preserves the long-range correlations. 

The expression 0 can be calculated for strong 
scattering in powers of the expansion parameter 
Ejj/p (El,: bandwidth), combined with a mean- 
field approximation as the starting point for the 
expansion. It has to be chosen such that the ex¬ 
pansion is convergent. 

Mean-field approximation: The (unnormalized) 
distribution density det C in Eq. 0 is approxi¬ 
mated by a constant phase arj ~ dj and the ex¬ 
act value is approached systematically by a conver¬ 
gent expansion in terms of the phase factor fluctu¬ 
ations — e“J. First, it should be noticed that 


with TTj = Tr2{(7jhphp^Q) for j = 0,1,2, 

hp = ao 2ip{Hp - ipao)~^ , (12) 


and with the 2D unit vector s = (cos A, sin A). 
Then Eq. 0 reads in this mean-field approxima¬ 
tion 


Prr' 





E 


CpiAi) 


(13) 


The pole of l/C'p(A;) with respect to e gives, af¬ 
ter the analytic continuation e ^ iui, the effective 
dispersion ojp of the new mode. This pole depends 
strongly on the details of the Hamiltonian compo¬ 
nents with Pauli matrices cti and cr 2 . It is impor¬ 
tant that we are only interested in the long range 
regime of Prr'fie)', i.e., in the behavior for small 
momenta. Therefore, we focus on hp for small p, 
and according to the expressions in m and 
this means that we need s ■ Pip for small p. In the 
following we calculate m for three different mod¬ 
els, namely 2D massive Dirac particles, 3D Weyl 
particles and particles on a square lattice with tt 
flux, whose low energy Hamiltonians behave like 
p[ a ■ p for small p. 

2D Dirac particles: For massive 2D Dirac 
fermions on a torus the Hamiltonian reads E[p = 
a ■ p a^m (0 < pj < 2 tt). Then we obtain 

~ 2ip . _x 

= CTo + ^5—- „ , (cr • p + cram -I- laap) 

and s = (cos(A-1-0), sin(A-1-0)) for Cp{A), which 
depends on the angle 9 = ej:ctan{2mp/(m? — ff')). 
For e ^ 0 we obtain 

1 _ (p^ + nP + ff)/2p 

Cp{A) 2e + is ■ p + pp"^/+ p'^) 
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random scattering. 


whose pole give the dispersion of a ray mode 

Wp \s-p-iDp^ with D = ^ . (15) 

The second term of Wp is imaginary and de¬ 
scribes damping. The mean-field condition 
maxA /plog(|C'p(A)|) is solved by A -I- d = 
0,7r/2, TT, 37r/2. Thus, four rays are created by 
strong random scattering with unit vectors si ,3 = 
(±1,0) and S 2,4 = (0, ±1) (cf. Fig. [T]). The lin¬ 
ear dispersion without damping is shown on the 
right-hand side of Fig. [51 

The dispersion (fT51) can also be understood as a 
Fokker-Planck dynamics with drift velocity s [l^ . 
This can be seen when we calculate the transi¬ 
tion probability from the expression (11411 within 
a Fourier representation in time and space. By as¬ 
suming that the propagation starts from the initial 
site 0 we get for the transition probability in Eq. 
@ for each of the four ray modes 


Pp{t) 




1 


(r ± TSif 


(16) 


with the effective diffusion coefficient D, which is 
the damping coefficient in Eq. (USD. This result de¬ 
scribes diffusion with a propagating (drifting) cen¬ 
ter in the direction of s). The damping term in 
the dispersion (USD leads to an additional diffusion 
away from these propagating centers. Diffusion 
slows down as the scattering rate ry is increased. 
This indicates that strong random scattering con¬ 
fines (or localizes) the movement of the particle to 
rays along the unit vectors s). 

3D Weyl particles: The Hamiltonian for 3D 
Weyl particles with 3D momentum p reads Hp = 
p- a. This implies for the propagator 


FIG. 2: Linear Dirac dispersion (left) and the linear 
dispersion of the four ray modes (right). The latter 
consists of four intersecting planes whose slopes agree 
up to a factor 1/2 with those of the Dirac dispersion. 


Then the dispersion of the ray modes reads 

]j^s ■ p - Wp^ . (17) 

Since the mass is zero in this case, the damping 
coefficient is D = l/2?7. Although the Hamiltonian 
is isotropic in 3D, there are again only four ray 
modes on the x-y torus. 

Square lattice with tt flux: The square lattice 
is bipartite. Therefore, the tight-binding Hamil¬ 
tonian for particles on a square lattice, which de¬ 
scribes nearest-neighbor hopping matrix elements, 
can be written in terms of two sublattices with co¬ 
ordinates (rj) (j = 1,2 and r are the coordinates 
on the sublattice with j = 1). Including the Peierls 
phase factors for the tt flux with ±i in x-direction 
and ±1 in y-direction, as indicated in Fig. [31 the 
corresponding Hamiltonian reads 

Hp = a ■ P + a^m, Pj = t sin akj . (18) 

This is a discretized 2D Dirac Hamiltonian with a 
lattice constant a and the hopping energy t. The 
flux through a plaquette of the square lattice is — 1, 
which is related to a phase tt. There is also a term 
with mass m that represents a staggered potential 
with m (—to) on sublattice j = 1 (sublattice j = 
2). This Hamiltonian also obeys property ([TD as 
described for the 2D Dirac Hamiltonian. The form 
of this Hamiltonian allows us to apply directly the 
results of the 2D Dirac Hamiltonian that gives 


/^ -. . -N 

CTo ± n , -2 (CT-P + icroV) 
p^ + rj^ 


from which we obtain 


1 _ ±? 7 ^)/ 2?7 

(7p(A) 2e + is-p + p^/fi ' 


1 {P"^ + m? + rf) / 27] 

C'p(A) 2e ± is • P ± fjP'^Km'^ + ff) 
This implies the dispersion 

ojp ^ —s ■ P — iDP^ 

^ 2 


( 19 ) 
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FIG. 3: Square lattice with tt flux as deflned by the 
tight-binding Hamiltonian na. The dots represent the 
sublattice with index j = 1 and the complex numbers 
are the Peierls phase factors for hopping from j = 1 to 
i = 2. 


for the ray mode on a square lattice with tt flux 
with the same directions as in the two previous 
cases. Here it should be noticed that there are 
small values for Pj near the four spectral nodes 
with akj = 0, ±7r. 

Discussion: The effect of strong random scatter¬ 
ing in models with a Weyl-Dirac Hamiltonian and 
for the tight-binding Hamiltonian on the square 
lattice with tt flux does not lead to conventional 
AL but creates four orthogonal ray modes. For pe¬ 
riodic boundary conditions the ray modes are par¬ 
allel to a torus or follow its circumference (cf. Fig. 
[T|) . The spontaneous creation of the four ray modes 
can be understood as an azimuthal localization at 
four discrete angles. In contrast to isotropic diffu¬ 
sion or isotropic Anderson localization, the states 
are confined to four orthogonal directions. This 
is similar to the uni-directional edge states in the 
quantum Hall effect or in photonic crystals with 
Faraday effect [^. The mechanisms for their cre¬ 
ation are rather different, though. Whereas the 
edge states are created in gapped photonic crys¬ 
tals at interfaces between regions with different 
Chern numbers Q, strong randomness creates the 
ray modes spontaneously in an isotropic system 
with strong random scattering. This makes them 
more accessible in real materials, provided that a 
band structure with linear dispersion exists for in¬ 


terband scattering. 

The result of the calculation can be summarized 
as a mapping of the momentum part of the Dirac- 
Weyl Hamiltonian to the ray-mode dispersion as 

a-p —?> -s-p — iDp^. (20) 

Thus, in all three cases we have considered here, 
the Pauli matrix vector a is replaced by the vector 
(l/2)s, whose components are scalars. Moreover, 
there is an additional damping term, whose coef¬ 
ficient D in Eq. (ESI decreases with an increasing 
scattering rate and increasing gap. The mapping 
(uni is reminiscent of the renormalization of the av¬ 
erage one-particles Green’s function as found, for 
instance, in the form of a complex self-energy in the 
self-consistent Born approximation. In the case of 
2D Weyl particles we get from random scattering 
the mapping 

a-p —>■ a ■ p + iaop . 

In contrast to the constant scattering rate the 
damping of the ray modes vanishes quadratically 
for small momenta. 

The fact that, for instance, photons can escape 
from a cloud of random scatterers is similar to the 
phenomenon of Klein tunnelling in a system with 
a potential barrier [H, [T^ . As in the latter case, it 
is crucial that there are two spectral bands. In our 
calculation this is reflected by the term s- tt in (fTTl) . 
It is sufficient for the generation of ray modes that 
this term vanishes linearly with the momentum. 

Tight-binding models on the square or the hon¬ 
eycomb lattice without flux do not develop ray 
modes because s • if (or equivalently s- H) is not 
linear for small momenta, as a calculation of the 
inverse propagator m indicates. However, the 
linear behavior oi s ■ H at small momenta might 
be obtained by breaking the time-reversal invari¬ 
ance: either by a periodic magnetic flux [l^ , by an 
additional spin texture on the honeycomb lattice 
M, or by the Faraday effect of a magneto-optic 
medium in a photonic crystal 
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